o : 
o ■ 

(N ■ 
<D ■ 

m : 
in : 



5-H 

CD 



C3 



I 

-a 
a 
o 

o 



> 
o 

oo 
m 

o 
l> 
o 



X 



Effects of an oscillating field on pattern formation in a ferromagnetic thin film: 
Analysis of patterns traveling at a low velocity 

Kazue Kudd3 and Katsuhiro Nakamura 
Department of Applied Physics, Graduate School of Engineering, Osaka City University, Osaka 558-8585, Japan 

(Dated: February 1, 2008) 

Magnetic domain patterns under an oscillating field are studied theoretically by using a simple 
Ising-like model. We propose two ways to investigate the effects of the oscillating field. The first 
one leads to a model in which rapidly oscillating terms are averaged out and the model can explain 
the existence of the maximum amplitude of the field for the appearance of patterns. The second 
one leads to a model that includes the delay of the response to the field and the model suggests 
the existence of a traveling pattern which moves very slowly compared with the time scale of the 
driving field. 

PACS numbers: 89.75.Kd, 75.70.Kw, 47.20.Lz, 75.10.Hk 



I. INTRODUCTION 

Rapidly driven systems have received considerable attention these days. Under a rapidly oscillating field, a state 
which is unstable in the absence of the oscillating field can be stabilized. One of the most simple and well-known 
examples is Kapitza's inverted pendulum and Landau and Lifshitz generalized the problem ll|. Their method has 
recently been applied to classical and quantum dynamics in periodically driven systems 0, Hf. The method is also 
applied to stabilization of a matter-wave soliton in two-dimensional Bose-Einstein condensates without an external 
trap 

Magnetic domain patterns in a uniaxial ferromagnetic thin film, which usually show a labyrinth structure, exhibit 
various kinds of structures under an oscillating field. For example, the labyrinth structure changes into a parallel-stripe 
structure for a certain field @, [|[. In some other cases, several types of lattice structures can appear [{|. 

In this paper, we develop effective theories for slow motion of magnetic domain patterns under a rapidly oscillating 
field. Especially, we focus on traveling patterns as an example of slowly moving patterns. So far, there were few 
effective theories to describe such a slowly traveling pattern under a rapidly oscillating field. In experiments on a 
garnet thin film, we can observe a parallel-stripe pattern traveling very slowly compared with the time scale of the 
field in some cases [lb] ]. A traveling mazelike pattern like Fig. [T]is also found in our numerical simulations. 

Although traveling patterns appear in various kinds of systems, most works about them have been limited to 
the systems in the absence of an oscillating field [HI, E2, EH, E3, EE SB, E3, EH ■ The mechanism of such traveling 
patterns in one-dimensional (ID) systems was intensively studied in the 1980s as drift instabilities or parity-breaking 
instabilities [1]], EH, EE EH- In Ref . [Till . secondary instabilities were discussed for several similar equations. By 
contrast, the authors of Refs. [H, EE Ej] gave no particular equation at first, but they considered symmetries of the 
system and assumed the form of the solution before deriving their equations. Almost 10 years after those papers, Price 
studied traveling patterns in 2D scalar nonlinear neural fields where the nonlinearity is p urely cubic and discussed 
constraints on the neural field structure and parameters to support traveling patterns [l5[ • He suggested that Swift- 




FIG. 1: Snapshots of a traveling pattern under an oscillating field after (a) 5000 cycles, (b) 10000 cycles, and (c) 15000 cycles. 
The whole pattern is traveling to the left. The details are described in Appendix lAl 
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Hohenberg-type models will not support traveling patterns. However, moving patterns were actually observed in 
numerical simulations for a general complex Swift-Hohcnbcrg equation in Ref. [la ] . Our model has similar properties 
to those of the (real) Swift-Hohenberg equation. We will not use the method of Ref. [l5[ but that of Ref. to 
explore the existence of a traveling pattern in a ferromagnetic thin film. 

In fact, recently, domain walls under a rapidly oscillating field have been studied [HI, [2(|. Michaelis et al. discussed 
the effects of rapid periodic oscillation of parameters in a Ginzburg-Landau (GL) equation by applying a multiscale 
technique and derived the averaged GL equation (l9j ; Kirakosyan et al. derived the averaged Landau-Lifshitz equation 
by employing the multi-time-scale expansion technique [201 ]. In those papers, they took into account higher harmonic 
oscillations. Although their methods cannot be directly applied to our model, our methods correspond to the lowest 
orders of their multi-time-scale expansions. 

Our model is a simple 2D Ising-like model (see Refs. [U, [23|, [24[ , and references therein), which has been used 
to simulate magnetic domain patterns. The numerical results simulated by the model show very similar properties to 
experimental ones [2l|,[22,[23|,l24|]. We consider a scalar field 4>(r), where r — (x, y). The positive and negative values 
of 4>(r) correspond to the up and down spins, respectively. The Hamiltonian consists of four energy terms: Uniaxial- 
anisotropy energy ff a m, exchange interactions ffj, dipolar interactions ffdi, and interactions with the external field 
-Hex- The anisotropy energy is given by 

fU.M.M. (1) 



It implies that the anisotropy energy prefers the values 4>(r) = ±1. The exchange and dipolar interactions are 
described by 

Hj = [Jl*pi (2) 



and 

ffdi = 7 / drdr'<f>(r)</>(r')G(r,r'), (3) 



respectively. Here, G(r,r') ~ \r — r'|~ 3 at long distances. The exchange interactions imply that <f>(r) tends to have 
the same value as neighbors. On the other hand, the dipolar interactions imply that <fr(r) tends to have the opposite 
sign to the values in a region at a long distance. Namely, ff j and may be interpreted as short-range attractive and 
long-range repulsive interactions, respectively. Their competition leads to a domain structure with a characteristic 
length. The term from the interactions with the external field is given by 

H ex = -h(t) fdr<p(r). (4) 

Here, we consider a spatially homogeneous and rapidly oscillating field: 

h(t) — hosinut. (5) 

From Eqs. dU) — d3J) , the dynamical equation of the model is described by 

d<P(r) SjH^ + Hj + ff di + ffex) 

dt ° 5(f>(r) 

= Lo(a[0(r)-0(r) 3 ]+/3V 2 0(r)- 7 / dr' <j>{r')G{r, r') + h(t)\. (6) 



Hereafter, we fix Lq = 1 and give the parameters a, /?, 7, and ho as positive values. 

In this paper, we propose two approximation methods to obtain the dynamical equation for slow motion. In both 
methods, we apply a part of Kapitza's idea that the dynamics under a rapidly oscillating field can be separated into 
a rapidly oscillating part and a slowly varying part jlj. In Sec. [Til we derive the model whose rapidly oscillating part 
is averaged out on the basis of the Kapitza's idea about the time average of the fast motion. In Sec. Illl[ we derive 
another model for the slow motion, considering the delay of the response to the field instead of taking a time average. 
After the derivation of the models, the instabilities of traveling patterns are investigated in both Sees. UTI and IIIII We 
discuss the details about the existence of a traveling pattern in Sec. IIVI Conclusions are given in Sec. [V] 
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II. TIME-AVERAGED MODEL 



First of all, we assume that the variable <fi(r) can be separated into two parts: 

<P(r,t) =$(r,t) + «£„(*). (7) 

Here, $(r,i) is a slowly varying term and 4>o(t) is a rapidly oscillating space-independent term. Substituting Eq. ([7)1 
into Eq. ©, we obtain 

+ = « [(#(r) + fo) - (*(r) + O ) 3 ] + /3V 2 ($(r) + fo) - 7 / dr'($(r') + o )G(r, r') + fc(t). (8) 

Let us consider only the rapidly oscillating space-independent part; then we have 

4> = a{4>o - 4>l) - 70o / dr'G(r'> 0) + h(t). (9) 

Here, we define G(r,0) = l/|r| 3 . Then, the integral in Eq. ([9]) is a constant, cto: 

r°° dr 

a Q = 2tt / — . (10) 



Here, d is the cutoff length to prevent the divergence for d —> 0. It is also interpreted as the lower limit of the dipolar 
interactions. The solution of Eq. ([9]) should have the following form: 

0o = posm(ujt + 6), (11) 

where S is a phase shift which comes from the delay of the response to the field. Substituting Eq. pT|) into Eq. © 
and omitting high-order harmonics (i.e. sm3cjt), we have 

3 

ujpo cos(ujt + S) = rjoPo s'm(ujt + S) — -^^Pq sin(uit + 5) + ho sinwi, (12) 
where rjQ = a — jao. From Eq. (|12[) . a pair of simultaneous equations is obtained: 

— ujpasmS = ^i]q - ^ a P^j Pocos5 + h , (13a) 
upoCosS = ^Tj — ^ap^J po sin<5. (13b) 

Eliminating 6 from Eq. (fTB"]) . we obtain a cubic equation of X = p\: 

l a 2 x 3 _ 3 a%x2 + (w a + %2)x = ^ (14) 

Therefore, po can be evaluated from Eq. (|14|) if the parameters a, ry , w, and /io are given. 

Now, let us think about the slowly varying part. After substituting Eq. (fTTj) into Eq. ([8|), we average out the rapid 
oscillation, closely following Kapitza's idea. Then, we obtain an equation for slowly varying domain patterns: 



0$(r) 
dt 



a ($(r) - $(r) 3 ) - ^ap 2 $(r) + /3V 2 $(r) - 7 / dr'$(r')G(r, r'). (15) 



The second term on the right hand side of Eq. (fT5j) is an extra term due to the time average. This term is essential 
to explore the effects of the rapidly oscillating field. 

On the basis of Eq. I|15p . we will analyze the possibility of patterns traveling at a low velocity. Let us first choose 
the most simple moving stripe- type solution for Eq. () 1 5[) : 



$(r, t) = A (t) + Ax(t) sin(kx + b(t)). (16) 

Substituting Eq. (|16|) into Eq. 1|15|) and omitting high-order harmonics, we have 

Ao + A\ sm(kx + b) +bAi cos(kx + b) = rj' Q Ao + r)[Ai s'm(kx + b) 

3 3 

Aq + 3AlA x sin(kx + b) + ^A A\ + -A\ sui{kx + b) . (17) 
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Here, 

Vo = f 1 - t;Po) a ~ 7«o, (18a) 

V'l = ^l-^p^a-/3fc 2 -7(a -aifc), (18b) 
with ao given by Eq. (fTU|) , a\ — 2tt, and k = \k\. Equation (jTTJ) leads to the following equations: 

A = ri' A - a (a 3 + ^A A^j , (19a) 

Ai = ii x A x - a \2>AlAi + -A? J , (19b) 

6 = 0. (19c) 
Equation (j!9c|) implies that the phase b(t) in Eq. IT6|) shows no time dependence and that there is no traveling pattern 



with the simplest form like Eq. |JJ 

Next, let us consider a more generalized solution by incorporating the second harmonics: 

<Z>(r,t) = A (t) + A 1 (t) sm(kx + b{t)) + A 21 cos[2(kx + b(t))] + A 22 sm[2{kx + b(t))}. (20) 

Substituting Eq. (fSU|) into Eq. (|15[) leads to the following equations: 

A = v'qAo -a\A% + ~A Al + ^A A 2 21 + ^A A 2 2 - jA 2 A 21 ^j , (21a) 

Ax = r&A! - a (jA* + 3AqAi + ~AxA 2 21 + \a x A\ 2 - 3A A X A 21 ^ , (21b) 

/3 3 3 3 \ 

Mi = r( 2 A 21 - a i-A 3 21 - -A A\ + 3A 2 A 21 + -A\A 21 + -A 21 A 2 22 - 6A a A 2 22 j , (21c) 

/3 3 3 \ 

A 22 = r( 2 A 22 - a f -A 3 22 + 3A 2 A 22 + -A 2 21 A 22 + -A 2 A 22 + 6A A 21 A 22 j , (21d) 



and 



Here, r]' and rf x are given by Eq. (fTS)) . and 



b=-3aA A 22 . (22) 



rfi =11- -Pi) a - Af3k 2 - 7 (a - 2a ± k). (23) 

This time, Eq (|22|) implies that there can be a traveling pattern if A$ ^ and A 22 ^ 0. 

Now let us find a stationary point (SP) of Eq. (|21[) where Ao = or A22 = 0, and examine its linear stability. If 
the SP is unstable and both A and A 22 grow from zero, the pattern can start to travel. For the parameter values 
used to obtain Fig. [TJ however, there are no SPs except for ones with Aq = A 2 \ = A 22 = 0. We should note Ax = 
or A\ = 4^ /3a at the SPs with Aq = A 2 \ = A 22 = 0. Since A\ must be real, rj{ > 0. Namely, 

Po < -5— [a — fik 2 — j(a — aifc)] . (24) 

This condition gives an estimate of the maximum value of the field amplitude ho to observe a nonuniform pattern, 
as ho proves to be a monotonic function of po f° r the parameter values in Fig. [TJ In other words, if Eq. (|24[) is 
not satisfied, the only SP is (Aq , A\ , A 2X , A 22 ) — (0,0,0,0), which means that no pattern appears. At the SPs with 
^0 = A 2 i = A 22 = and A 2 = irji/Sa, the Jacobian of Eq. (f2T|) becomes 

/ Vo - ^'1 i rfy \ 

3 = H 11 V ' 2 - 2 V [ ' (25) 

V V2-H/ 

The real parts of eigenvalues of Eq. ([25]) are Ai = — 2rj[, A 2 = A 3 = ^(r)' — 4r][ + r) 2 ), and A 4 = rj 2 — 2r][. Note that 
A x is always negative. The others (A 2 , A3, A4) also prove to be negative when k ~ 1. In fact, the most preferable 
wave number of domain patterns is k = 1 for the parameter values in Fig. [JJ (see Ref. [24J for details). Therefore, the 
present SPs are stable and we cannot expect a traveling pattern in this case. 
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III. PHASE-SHIFTED MODEL 



In this section, we consider another equation for slowly varying domain patterns instead of Eq. (|15[) . We begin 
with Eqs. (|7))- (fH|) again, but we will not take a time average. Instead, we take the delay of the response to the field 
into consideration. Substituting Eq. (fTTj) into Eq. ([5]), we consider the equation as a discrete-time equation which is 
valid at t = (2ir/uj)n with integers n. Then, we regard the discrete time as continuous. This procedure is justified 
when the field oscillation is rapid enough compared with the time scale of the slowly varying part. It is as if we take 
a sequence of snapshots at t = (2ir/uj)n and take it as a movie. In fact, our numerical results in Fig. [T]are obtained 
by taking these kinds of snapshots. We thus obtain a new equation for slowly varying domain patterns: 



<9$(r) 
dt 



a(l - 3pl sin 2 <5)$(r) + /3V 2 $(r) - 7 / dr'$(r')G(r, r') - a$(r) 2 ($(r) + 3p sin 5) + C, 



where 



C = ryoPo sin 8 — ap^ sin 3 8 — uipo cos S, 



(26) 



(27) 



with po and 8 evaluated from Eq. (fT3| . Equation (|26| has two extra terms due to the phase shift 8 except for the 
constant C. One is linear and the other is nonlinear in The extra nonlinear term has an important role in discussion 
of the existence of a traveling pattern. 

Now, let us consider the stability of a traveling pattern on the basis of Eq. (|26|) . When the simplest form, Eq. (|16l) , 
is substituted into Eq. (|26|) . we obtain the same result as Eq. (|19c[) . Therefore, we proceed to choose the extended 
solution, Eq. ([20|) . Substituting Eq. ([20]) into Eq. ([26]) leads to the following equations: 



and 



A = ij a A + C 



-a 



A 2 (A Q + 3p sin 8) + -(A + Po sin 5){A\ + A 2 21 + A 2 22 ) - -A\A 21 



A x = fjiM 



-A 



Ah) 



3A A 1 {A + 2p sin 5) - 3A 1 A 21 (A + p sin 6) 



A21 = f] 2 A 



21 



-A^ + - A 21 (2Al + Al 2 ) + 3 A A 21 {A + 2p sin 6) - -{A\ + AA 2 22 )(A Q + p sin 5) 



V2A2 



3 A* 



-A 22 {2A\ + A 2 21 ) + 3^0^22(^0 + 2p sin 8) + 6^21^22(^0 



posm.8) 



(28a) 



(28b) 



(28c) 



(28d) 



Here, 



m 
m 



(1 
(i 
(i 



b = -3a(A + po sin<5)yl22- 

3pp sin 2 8)a — 700, 

3pp sin 2 8)a — f3k 2 — 7(00 - 

3po sin 2 8)a - 4/3fc 2 - j(a 



a\k), 
- 2aifc). 



(28e) 

(29a) 
(29b) 
(29c) 



Equation (|28c[) suggests that there can be a traveling pattern if both Aq + po sin 8 ^ and A 22 ^ are satisfied. 

Now let us think about the SPs of Eq. ([28]) where Ao + po sin 6 = or A22 = 0. For the cases with k ~ 1 and the 
parameter set used in Fig. [TJ we find that there are no SPs with Ao + pg sin 8 = 0. Therefore, we concentrate on SPs 
with A22 = 0, where Eq. (|28|) leads to the following equations: 



= fjoAo + C -a 
= fjiAi - 
= 772^21 - a 



A 2 (A + 3p sin 8) + - (A + Po sin 8) (A\ 



A 2 21 ) 



A\A 21 



3^ Ai(A + 2p Q sin 8) - 3A 1 A 21 (A a + p a sin 8) 

3 



3 3 < 

-A 3 21 + -A 2 A 21 + 3A A 21 {Ao + 2p n sin 8) ~ -A 2 {A + p sin 8) 



(30a) 
(30b) 
(30c) 
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Here, we note that A\ should not be zero. When A\ ^ 0, Eq. (|30b|) leads to 



A\ = -^--2[Al 1 + 2 A a (A + 2 Po sin 5) - 2A 21 ( A + Po sin 5)] 



(31) 



Substituting Eq. (|3ip into Eqs. (|30ap and (|30c|) . we obtain a pair of nonlinear simultaneous equations for Aq and A21, 
which can be solved numerically. 

At those SPs, the Jacobian of Eq. ([28| is 



J = 



I Jn J12 J13 

2J\2 J22 J23 

2^13 J23 J33 

V J 44 



(32) 



with the elements given by 
J11 
J12 

Jl3 

J22 
J23 
J33 



3A 2 + 6A oPo sin 5 + - {A\ + A 2 2l ) 



770 - a 



-3aAi ( A + po sin S - - A 2 i 



-3a 



(Aq + Po sm5)A 21 - -A( 
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-A\ + -A\ x + 3A {A + 2 Po sin 6) - 3A 21 (A + Po sin S) 



J44 = m - a 



rji-a 
3aAi (A + po sin S - A 2 i 
fj 2 -a 

"3 
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4^21 + 2 A i + 3 MAo + 2p sin 5) 



-^(2A\ + A 2 21 ) + 3A a (A + 2p a sin 5) + 6A 21 {A a + p sin 6) 



(33a) 

(33b) 

(33c) 

(33d) 
(33e) 
(33f) 

(33g) 



Equation (|32|) is a block-diagonal matrix. We can evaluate the real parts of the eigenvalues, Ai, A2, A3, for the 
upper-left 3x3 matrix as well as A4 = J44. The dependence of the real parts of the eigenvalues on the field amplitude 
h is shown in Fig. O Here, we take the values of the parameters, a, /3, and 7, used in Fig. [TJ For k = 1.0, all the real 
parts of the eigenvalues (Ai, . . . , A4) are always negative. In other words, the SPs are stable and a traveling pattern 
cannot appear. For k — 0.83, however, extra SPs appear in the region between ho — 0.35 and ho ~ 1.5. In that 
region, A4 and one of the other three (A 1 ,A2,A 3 ) are positive. Incidentally, it is confirmed that A Q + posing 7^ in 
the region. This result suggests that a traveling pattern can appear in a certain region of the field when k — 0.83. 

Using the above analysis, we show a stability diagram in Fig. [3] In the unstable area, where there is a branch with 
positive A4, a traveling pattern can appear. In the stable area, where all the branches of A4 have negative values, 
it cannot appear. The values of the parameters a, (3, and 7 are the same as ones in Fig. [TJ The characteristic wave 
number fco depends on the ratio of (3 and 7 (see Ref. 24| for details), and fco = 1 in our case. Though it is expected 
that k ~ ko, the actual characteristic length in the simulations is larger than 27r/fco- I n other words, fc < fco in the 
actual numerical results, although a domain pattern with a small fc is not always realistic. Incidentally, the field larger 
than ho — 1.5 may be meaningless since domain patterns should vanish under a strong field. 



IV. DISCUSSION 



The results in Fig. [5] suggest that a traveling pattern can exist for fc < 0.83 but not for fc = 1.0. This can be 
interpreted as meaning that a traveling pattern should be a little fat. In fact, the actual wave numbers of domain 
patterns in our numerical simulations are a little less than fc = 1, although fc = 0.83 seems too small. We can say 
that this fact partly supports the theoretical results given here. 

We have used very simple approximations, i.e. perfect parallel-stripe structures without any distortion, to investigate 
the instabilities of a traveling pattern. That may be one of the reasons why the present analysis has suggested a 
traveling pattern with fc smaller than that of the numerical results. In our simulations, the traveling patterns do not 
have a perfect parallel-stripe structure. If more complex and better approximations are employed, the actual traveling 
patterns exhibited by numerical simulations may be better explained. 



7 



< 

CO 

< 
< 
< 



(a) 



Ai, A2, A3 x 
A4 • 



xxxxxxxxx xxxxxxxxxxx 



xxxxxxxxxxxxx x x x 



— m« »«»00»00OO<XXXXXX X xxx xxxxxxxxxxx 



0.5 



ho 



1.5 



< 
< 
< 
< 



(b) 



^^^SCKXXXXXXXXXXXX X X X X x 

xxxxxxxxx^ 



xxxxxxxxx X XXXXXXXXX X 



^^ 1 **««««<XXXXXXXXXXX X X X x x x 

smW8S88xxx»«><xxxxxxxxxxxx x x'x'fSt******* 



xxxxxxxxx x xxxxxxxxx X 



0.5 



1.5 



FIG. 2: The dependence of the real parts of the eigenvalues of Eq. (I32[l on the field amplitude ho: (a) k — 1.0 and (b) k = 0.83. 



1.5 



ho 



0.5 



00000 



0000 
0000 



OOOOOOOOOOOOOXXXXXXXXXXXXXXXXXX 



00000 
00000 



oooooxxxx 
oooooxxxx 



xxxxxxxxxxxx 
xxxxxxxxxxxx 



00000000000000 
0000 

0000 
0000 



X X X X X X 



X X X X X X 



unstable 



oooxxxxx 



stable 



X X X X X X X 



ooooooooooooooxxxxxxxxxxxxxxxx 



00000 
00000 
00000 
00000 
00000 
00000 



00000 
_ 00000.. 
oooooooxx:. 

oooooooxxxxxxxxxxxxxxxxx 
oooooooxxxxxxxxxxxxxxxxx 



X X X X X X 



X X X X X X X 



X o o o o 

x x Q o o „ 

X X x o o o o 

X X X v o o o 

x it x 9 q O O 

8 8 g 



o 
o o 
o o 



00000 
00000 
00000 
00000 



O O O O O X X 



X X X X X X 



xxxxxxxxx 

■ X X X X X X 



OOOOOXXXXXXXXXXXXXXXXX 



X X X O 



X X X X X 



o o 

o o 

o o 

o o 



00000 

O O O O v 

0000: 



O O O O X X X 



xxxxxxxx 



X X X o o 



O X X X X X X 
X X X X X 



X X X X X X 



X X X X X X 
X X X X X X 



0.7 



g 3 3 $ 
0.75 







0.8 



3 3 3 3 
0.85 

k 







3 3 3 
0.9 



3 3 3 
0.95 



3 3 3 3 

1 



FIG. 3: Stability diagram in k-ho space. In the unstable area (red-circle points), a traveling pattern can appear. In the stable 
area (green-cross points), the pattern cannot travel. 



In experiments, the perfect parallel-stripe structure is a realistic pattern. However, as mentioned above, the 
condition for a traveling pattern is tight even for such a simple structure. Traveling patterns with a more complex 
structure can be observed in experiments and the conditions of their appearance would be more complex than the 
present case. In any case, it is sure that a traveling pattern cannot appear without a rapid oscillating field. 
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V. CONCLUSIONS 

We have proposed two ways to describe magnetic domain patterns moving slowly under a rapidly oscillating field. 
One gives a model in which rapidly oscillating terms are averaged out. The time-averaged model can explain the 
existence of the maximum values of the field where non-uniform domain patterns are preferable. The other gives a 
model which includes a phase shift as the delay of the response to the field. The phase-shifted model suggests the 
existence of a traveling pattern which moves very slowly compared with the time scale of the field. These two models 
have both merits and demerits. In other words, the approximations to be employed depend on the phenomenon under 
consideration. We should choose a method suitable for the analysis of the phenomenon to be investigated. 

Although we have focused on a traveling pattern in this paper, these two methods are promising for applying to 
many other domain patterns under a rapidly oscillating field. 
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APPENDIX A: NUMERICAL SIMULATIONS 

The numerical procedures for time evolution are almost the same as those of Refs. [23L |24|. For time evolution, we 
use a semi-implicit method: The exact solutions and the second order Runge-Kutta method are used for the linear and 
nonlinear terms, respectively. For a better spatial resolution, a pseudo-spectral method is applied. In other words, 
we numerically calculate the time evolutions of Eq. ([6]) in Fourier space: 

^ = c# - 3 ] fc - {(3k 2 + jG k )^ k + h(t)S k , (Al) 
at 

where denotes the convolution sum and G k is the Fourier transform of G(r,0). Since we defined G(r,0) = l/|r| 3 , 

Gk = a-o - aik, (A2) 

where k = |fc| and 

r°° dr 

a = 2-K — , 01 = 2tt. (A3) 
Jd r 2 

In the simulations, we set d = tt/2, which results in ao = 4. 

In Fig. [l] the parameters are given as a = 2.0, f3 = 2.0, and 7 = 2/3/ai = 2/ir. The frequency and amplitude of 
the field are u> — 2ir x 5 x 10~ 2 and ho = 0.8, respectively. The simulations are performed on a 128 x 128 lattice with 
periodic boundary conditions. The snapshots in Fig. [1] are the domain patterns at (a) 5 x 10 3 T, (b) 10 x 10 3 T, and 
(c) 15 x 10 3 T, where T — 2-k/uo. If the amplitude is larger (for example, ho = 0.9, 0.95, etc.), we can see a traveling 
pattern with a different structure moves to a different direction. 
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